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Abstract 

We present and prove a theorem of matrix analysis, the Flavour Expansion Theo¬ 
rem (or FET), according to which, an analytic function of a Hermitian matrix can be 
expanded polynomially in terms of its off-diagonal elements with coefficients being the 
divided differences of the analytic function and arguments the diagonal elements of the 
Hermitian matrix. The theorem is applicable in case of flavour changing amplitudes. At 
one-loop level this procedure is particularly natural due to the observation that every loop 
function in the Passarino-Veltman basis can be recursively expressed in terms of divided 
differences. FET helps to algebraically translate an amplitude written in mass eigenba- 
sis into flavour mass insertions, without performing diagrammatic calculations in flavour 
basis. As a non-trivial application of FET up to a third order, we demonstrate its use in 
calculating strong bounds on the real parts of flavour changing mass insertions in the up- 
squark sector of the MSSM from neutron Electric Dipole Moment (nEDM) measurements, 
assuming that CP-violation arises only from the CKM matrix. 
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1 Introduction 


Within the general framework of a perturbative Quantum Field Theory (QFT), the standard 
strategy followed when calculating physical transition amplitudes, is to express the Lagrangian 
density in a particular field basis, commonly referred to as mass eigenstate basis. Contrary 
to any other possible choice, only in this case the states of the theory correspond to physical 
particles with definite mass and symmetry charges. Up to possible mass degeneracies, this 
basis is unique and is characterized by the absence of quadratic mixing terms between different 
mass eigenstates. Furthermore, all parameters of the Lagrangian in this basis, are physically 
observable, in the sense that all masses and couplings can in principle be obtained by a 
suitable experiment. After having set the Lagrangian to the mass eigenstates fields basis, one 
can then deploy the standard QFT machinery and set the Feynman rules in order to calculate 
transition amplitudes for any physical process. 

Nevertheless, in the vast majority of the models we are interested in, masses are typically 
generated or affected by a symmetry breaking mechanism. In this case another basis is 
physically meaningful as well. This is the basis where the Lagrangian exhibits explicitly 
the properties of the initial symmetry, and the states correspond to eigenstates of a larger 
symmetry group. We will refer to these eigenstates, for gauge bosons and collectively for 
fermions and scalars with family replication, as gauge and flavour eigenstates respectively, 
although in our definition for the latter there is no implicit requirement of an underlying 
flavour symmetry. In this sense the flavour eigenstate basis in many models can be considered 
in practice arbitrary, constrained only by the other symmetries of the initial Lagrangian, 
i.e., gauge symmetry, supersymmetry, etc. The transformation from the initial basis to the 
mass eigenstate basis, which is still the physical basis of the theory, is performed with mass 
diagonalization involving unitary transformations and field redefinitions, typically leaving a 
physical imprint on the parameters of the mass eigenstate theory. In the Standard Model 
(SM) [iJI^ this effect is displayed in the gauge sector through the weak mixing angle and 
in the fermion sector through the CKM and PMNS matrices of charged currents. 
However, even in this very successful model, the CKM or PMNS parameters along with the 
fermion mass eigenvalues are insufficient to determine unambiguously the flavour eigenstate 
basis. 

Although the mass eigenstate basis of a perturbative QFT is the natural basis for calcula¬ 
tions of physical processes, some effects typically related to the symmetries of the Lagrangian 
before symmetry breaking are better understood in flavour basisQ Therefore for a qualitative 
analysis of such effects, it is often useful to have our expressions in the latter basis0 The stan¬ 
dard strategy that has been employed up to date, is an approximate diagrammatic method 
commonly referred to as the Mass Insertion Approximation (MIA) |^[^. In this approach 
one defines the Feynman rules of the theory directly in flavour basis. The diagonal part 
of the flavour mass matrix is typically absorbed into the definition of (unphysical) massive 
propagators and the non-diagonal part commonly referred to as mass insertions is treated 
perturbatively, as part of the interaction Lagrangian which now possesses quadratic mixing 
terms. Due to the presence of quadratic interactions, besides the standard loop approxima- 

^Since in many cases the mass diagonalization of various sectors of the theory is independent of each other, 
one can also work in a mixed basis where some sectors are given in mass basis and others in flavour basis. In 
what follows, the basis we work can be easily identified from the context. 

^This is after all the basis that more naturally connects couplings and masses to high energies through their 
Renormalization Group Equations (RGEs). 
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tion of a perturbative QFT, there is an extra approximation for each diagram, appearing as 
an infinite series in terms of mass insertions, following the presence of a flavour propagator. 


In what follows, we present an algebraic treatment of transition amplitudes in mass eigen¬ 
state basis, leading directly to the corresponding amplitudes in flavonr basis, in the form of 
the MIA or of an equivalent expansion. In particnlar, we prove a theorem in matrix anal¬ 
ysis [^,[H , that we coin Flavour Expansion Theorem or simply FET, which says that an 


analytic function of a Hermitian matrix can be expanded polynomially in terms of its off- 
diagonal elements with coefficients being the divided difference of the analytic function and 
arguments the diagonal elements of the Hermitian matrix. At one-loop level, this expansion 
is naturally related to the remarkable recursive properties of next order Passarino-Veltman 
(PV) function [12] being the divided difference of the previous one. We then argue that 
FET connects mass and flavour field bases amplitudes. The first non-trivial order in the 
expansion [cf. eq. (3.11)], and applications in MSSM flavour physics, have been presented 
in refs. [14[|15] but a formal mathematical proof to all orders was unknown until now. FET 
is especially useful when is used to evaluate higher order expansion terms, it is technically 
easier, elegant and superior to often tedious, time-consuming and thus prone to mistakes 
calculations of the diagrammatic MIA. We support our claims with a novel example towards 
the end of the article. 


More specifically, the paper is organized as follows: In Section]^ we present a warming-up 
example of a scalar toy-model in order to illustrate the relation between the calcnlation of 
flavour transition amplitudes in mass and flavour bases. Then, in Section]^ we formulate a 
general algebraic theorem dealing with the expansion of an analytic function of a Hermitian 
matrix, and, discuss its applications to rewriting flavour amplitudes with scalar and vector 
particles, from mass to flavour eigenstates basis. We extend our discussion to the case of 
amplitudes involving fermions in Section]^ In Section]^ we illustrate the developed technique 
on a (potentially) physical example, expanding the dominant gluino-squark contribution to 
the neutron Electric Dipole Moment and showing the importance of higher order terms. We 
conclude our results in Section]^ Finally, the formal proof of theorem formulated in Section]^ 
is given in Appendix A[ while in [Appendix B we derive the convergence criterion for the mass 
insertion expansion of the one-loop integrals. 


2 A warming-up example: flavour calculation techniques 

To set up a simple framework to introduce the standard techniques of flavour physics cal¬ 
culations, we consider a scalar field toy model composed of iV-complex charged scalar fields 
$/, with family replication, and an extra neutral, real, scalar field rj. The (squared) mass 
matrix, M^, and the Yukawa coupling matrix, Y, of the flavour eigenstates <h/, are necessarily 
Hermitian but not aligned in general. The Lagrangian density, will have the form]^ 

>Cflavour = - Yjj . , (2.1) 

where dots denote additional scalar field interactions which are irrelevant for the discussion 
below. Using the unitary rotation, 

‘h/ = Uii (f>i , (2.2) 

^Sum over repeated indices is always assumed in the text, unless stated otherwise. 
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Figure 1: Scalar self-energy —i'Eji{p) in the mass eigenstate basis. 
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Figure 2: Scalar self-energy —i'Eji(p) in the flavour basis. 


where U satisfies the condition 


(2.3) 


(2.4) 


U = = diag(m7,..., m%) , 

one can express the Lagrangian in terms of mass eigenstates (pi 

/Imass = (9^4) {d^(pi) - mj(p\(pi + ^(<9^p) {d^p) - - pij p (p\ (pj + ... 

where the transformed scalar couplings are identified as pij = Ujj Yjj Ujj. 

First, let us consider the “flavour” changing one-loop One-Particle-Irreducible (IPI) self¬ 
energy diagram of the mass eigenstates fields (pi, shown in Fig. Then using interactions 
from Lagrangian in eq. (2.4) leads to 


-iT.ji{p) = 


d^k 


(27r)^ {kP — mj){{k — pY — m^) 


2\ Vh 


(47r) 


yjuBoip-1711,171^) ya , (2.5) 


where the loop function Bo (and also Co, Do ,... below) is a PV-function defined in Section 3.3 


Next, we consider the corresponding one-loop diagram in flavour basis of eq. (2.1), employ¬ 
ing the standard diagrammatic MIA approach. The massive flavour propagators are defined 
by absorbing the diagonal part of the flavour mass matrix M^, according to the decomposition 
into diagonal and non-diagonal parts, 


Mfj = {Mq)ii 5ij + M]j , Mfj = 0, (no sum over I) , (2.6) 

where d/j is the usual Kronecker symmetric tensor. The non-diagonal elements M‘jj, are iden¬ 
tified as the mass insertions, treated as perturbative couplings for the non-diagonal quadratic 
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interactions of flavour fields. The one-loop flavour changing self-energy of the flavour states, 
is thus represented as the infinite sum of the diagrams shown in Fig[^ By direct calcu¬ 
lation, and denoting = {Mq)kk, one obtains 


- i'^Jiip) = 


+ 


d^k YjkYli 


5kl 


+ 




(27r)^ {k — pY — \k'^ — M"^ {k? — M‘^){k‘^ — M^) 




(A:2 - M|.)(A:2 - M^Yk"^ - Ml) 




(2.7) 


(4^2) YjkYli x ^klBo{p', Ml^,mlj) + M^£Co(0,p; M|, m^) 


+ Ml^N^NLBoiO,0,P] Mj^, Mlf, MI,!!!^) -b 


( 2 . 8 ) 


which is essentially an expansion in terms of mass insertions. We should notice that although 
are not the squares of the physical masses, they are always real non-negative. This is a 
general property of the diagonal part of any semi-positive definite Hermitian matrix, including 
also any Hermitian (squared) mass matrix of a consistent QFT. 

As the external indices imply, the self energy diagrams in this example are not invariant 
under flavour rotations. One can formally uncover the explicit correspondence between the 
flavour and mass basis calculations by considering the flavour invariance of the time evolution 
operator, inside the corresponding S-scattering matrix element, with the relevant contractions 


= d'^x d'^y{p,J\^]{x)0^^l.y^(l)i{y)\p',I) , 


{p,J\^]{x) = 


^ipx 


Uj, , = , 

we derive the transformation rule for self-energies, 

Sj/(p) = Ujj Yjiip) 4 > 

which can be immediately generalized to the case of an arbitrary n-point amplitude. 


(2.9) 

( 2 . 10 ) 


Substituting in eq. (2.10) the explicit algebraic expressions for the self-energies we obtain 
an interesting result - the flavour rotation of the mass eigenstates loop-function is an expansion 
in terms of mass insertions in flavour basis (no sum over K, L), 

Uki Bo{p,mj,ml) = 6KLBo{p]MK,ml) + MKLCo{0,p;MK,Ml,ml) 

+ MI^n^nlBo{ 0,0,P', MI^, Mff, Ml,ml) -b ... . (2.11) 


This result, however, can be also obtained by a theorem of matrix analysis [cf. eq. (3.11)] 
stated in the next section, rendering diagrammatic calculations in flavour basis, similar to 
ones leading to eq. (2.8), obsolete. 


3 Flavour Expansion Theorem 


Eq. (2.11) has been obtained diagrammatically with the help of the transformation rule (2.10). 
In what follows we show that, such relations can be also obtained purely algebraically, allowing 
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for an easier transformation between mass and flavour basis calculations without the use of 
the diagrammatic MIA. In this section we formulate a relevant mathematical framework and 
a useful general theorem of matrix analysis. For brevity we refer to it as “Flavour Expansion 
Theorem” or just FET. 


Before formulating EET, it is worth noting that obtaining the relation (2.11) in a closed 
form without reverting to diagrammatic MIA expansion is not easy with the use of standard 
perturbation techniques. The simplest idea of expanding the mass eigenstates result in a 
Taylor series around some average mass uiq = ^ Y1k=i 


U^Kf{ml)U]K = UiK 


E 

_ n=0 


/H(, 


mn 


2\n 


n\ 


{rriK - mo) 


u: 


jK > 


(3.1) 


fails to reproduce correctly even the first non-trivial MIA term in eq. (2.11) - the higher terms 
of any order in Taylor expansion are proportional to a factor 


U,K{ml-ml)^U*i, = = [(m^ - , (3.2) 

which after expanding would explicitly contain terms linear in Mfj. 

Alternatively, a more consistent approach can be developed using the standard quantum 
mechanic perturbation theory, applied to mass matrix eigenstate problem, 

Mq + + ...) = {m‘^)f^ + X{m‘^)^P + ... + Ae^^^ + ...) ,(3.3) 


with A being the expansion parameter. By solving eq. (3.3) order by order we find mass 
eigenvalues mj and rotation matrices U = (ei,..., gn) as a series in A. Then the product, 
UiK can again be expanded in Taylor series in A parameter, with each term 

equivalent to the same order of MIA expansion. Such procedure easily restores the first terms 
in eq. (2.11) (see e.g. (^), but its combinatorial complication grows quickly and again it is 


hard to see how the higher order terms combine to form compact n-point loop functions, a 
situation which becomes even trickier in case of degenerate eigenvalues. 


3.1 Hermitian matrix function and Divided Differences 

We append below definitions that are relevant for presenting the Flavour Expansion Theorem. 

Definition 1. Let A be an n x n Hermitian matrix, diagonalized by a unitary matrix U to 
a real diagonal matrix D, through the transformation, 

UUU = D = diag(Ai,...,A„). (3.4) 

Also let f{x) be a real analytic function about zero, in a range / C IR, that can be expressed 
in terms of its Maclaurin series, as 

OO 

fix) = ^ CmX"" . (3.5) 

m=0 

Then, if all Xi E I, one can define a Hermitian matrix function f{A), as 

OO 

U/(D)Ut = ^ c^UD^Ut = ^c™A™ = /(A). (3.6) 

m=0 
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Definition 2. For any function f{x), one can define a set of multi-variable functions, f^^\ 
through a mathematical operation commonly referred to as divided difference. Divided differ¬ 
ences are defined recursively as 


/[°l(x) = /(x) , 


/W(xo,xi) 


f{xo) - f{xi) 

xo - xi 


(3.7a) 

(3.7b) 


/[^+^l(xo,. . .,Xk,Xk+l) 


/W(xo,...,,Xfc-l,XA;) -/W(xo,... ) Xk—lj Xk-{-l) 
Xk Xk-\-l 


(3.7c) 


It is easy to check that divided differences of any order k are always totally symmetric 
under the permutation of any set of respective arguments, Xj. Moreover, for an analytic 
generating functional f = f^^^ they also have a well defined degeneracy limit 


lim 


f^^\xo,...,Xk) 


1 (9™ 
m! 


/["-™l(e,x™+i...,Xfc), 


(3.8) 


applied to any set of respective arguments (m < k), as long as all arguments he in the domain 
of analyticity of f{x). 


3.2 The Flavour Expansion Theorem (FET) 


By making extensive use of the definitions above, we can formulate a general theorem, con¬ 
cerning a certain expansion of Hermitian matrix functions, in a form which can be directly 
applied to the unitary transformation of loop functions in flavour physics. 

Theorem 1. Let A be any nxn Hermitian matrix, decomposed as a sum of its diagonal and 
non-diagonal part, 

A = Ao + A, (3.9) 

where, by definition, 

(3-10) 

Aij = Aij, All = 0 , (/, J = 1,..., n) . 

Then, for any Hermitian matrix function /(A), satisfying the restrictions of Def.^ a given 
matrix element {I\f{A)\J) will be given by the expansion (no sum over I,J) 


f{Ad)ij = (5/j/(Ao) -h /t^'(Ao, Ag) Aij -h ^/[^'(Ag, Ag, A^^) Aik^Ak^j 


Ki 


AiKiAk.,^K2-^K2J A ■ ■ ■ ) 


(3.11) 


Ki,K2 


in terms of divided differences o//(Ao) [see Def.^ and the (non-diagonal) elements o/Aj^ 


Eq. 3.11 holds as long as the expansion in the RHS is convergent. The formal proof of this 
theorem, based on the notion of fully symmetrized polynomials and mathematical induction. 
Appendix A 


IS given m 


^Cases of degenerate eigenvalues and/or diagonal matrix elements are treated uniformly due to property 


(3.8 1 of the divided differences. 
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3.3 Divided Differences and Passarino-Veltman functions 


The natural connection between FET and the expansion of one-loop amplitudes is becoming 
striking when looking into the recursive properties of the loop functions in the Passarino- 
Veltman basis |12]. The general n-point one-loop functions can be defined as: 


(dvr) 




d^k 


... k^’- 


(27r)4 (A:2 - ml) nL2((^ + Pi H-h - ml) ’ 


(n > 2) 


(3.12) 


In the standard notation n = 2,3 ... functions are commonly denoted as B,C, .. .-functions. 

A useful property associates differences of integral functions of a certain order with integral 
functions of next order. In general case this relation has the following structure: 


PVn (Pi ■ ■■Pn-i-,m{ ... mi) - PV^ {pi. ..pn-i;m[^ ... mi) 


mt — m 


/2 


= PI4^i(0,pi.. . ..mi) 


/2 


(3.13) 






) - PV^{...pj-i...;...m)\..) 


/2 


2 / 2 
rrij — m- 


= PI4+i(...Pj-i,0...;...m^,m'^..), (j > 2) 

with X being any set of Lorentz indices of momenta in the numerator of loop integrandj^ 
Comparing eq. (3.13) with the definition (3.7) one can see immediately that the notion 
of divided differences is naturally implemented in the relations between multi-point one-loop 
integrals. Eq. (3.13) allow us to express the expansion (3.11) for one-loop amplitudes in a 
form in which it is obvious that it is not singular for degenerate diagonal matrix elements. 
Namely, every one-loop amplitude can be written as a linear combination of PV functions. 
Purthermore, each PV function can be expanded as, 


PV^^\...,A,...) = 5ijPV^^\...,Ai,...) + PV^^+^\...,Ai,Ai,...)Ajj 


J/j 


+ Aq , Aq ,.. .)AikAkj + 


(3.14) 


K 


where if necessary one should also specify momenta arguments as prescribed in eq. (3.13). 

Eor example, to make a connection between PET and the toy-model of the previous section 
we observe that if we make the following identifications 


D 


m 




fix) = /f°'(x) Bo{p,x,ml 


(3.15) 


we can immediately see that eq. (2.11) is a special case of eq. (3.11). In Section]^ we will 
present a highly non-trivial example application of the PET. 


^Additional discussion and more recursive relations for the various types of PV functions can be found in 
Appendix A of ref. 
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3.4 Applications and limitations of FET 


The FET formulated as a pure algebraic theorem can be directly applied to expanding a mass 
eigenstate result of any transition amplitude in any model involving particles associated with 
Hermitian mass matrices, that is scalars or vectors, even at higher loop orders. As we shall 
discuss in Section with some modifications, FET can be also applied to expanding the 
amplitudes involving fermions associated with non-Hermitian mass matrices. 

The purely algebraic expansion is usually significantly simpler than the more tedions and 
prone to mistakes diagrammatic MIA calcnlation, particularly in models with complicated 
flavour structure like MSSM. Another advantage of FET is that it can be easily implemented 
as an algorithm for symbolic manipulation programs, antomatizing the expansion procedure. 
However, the procedure has some limitations, particnlarly when is applied to such complicated 
functions as loop integrals. Three remarks concerning limitations of EET are summarized 
here: 


Remark 1. EET assumptions require the expanded amplitude to be analytic function of 
masses. This is not the case if external momenta are large and loop integrals may have 
branch cuts. For the external momenta in the vicinity of branch cnts a mass eigenstate 
calculation and use of numerical procedures is more appropriate. 


Remark 2. Flavour expansion in the r.h.s of eq. (3.11) may not converge or converge very 
slowly, in both cases again mass eigenstates basis and nse of numerical procedures is preferred. 


Remark 3. The UV-singularities do not pose a problem for FET. If they appear, they come 
from the loop integrals of positive mass dimension. The coefficients of poles of such integrals 
are dimensionless and flavour blind or proportional to positive powers of masses, so they can 
be evaluated in terms of flavour basis parameters without any expansion. 


Most practical applications of FET concern analyses of models of New Physics where 
loop particles are mnch heavier than the external states (being usnally the Standard Model 
fields). Thns, it is nsually sufficient to calculate relevant amplitudes in the approximation 
of vanishing external momenta, or by expanding the loop integrals in the external momenta 
and keeping only the first few terms of such an expansion. Since in these processes the loop 
integrals are real analytic functions of masses, branch cuts can never appear. Then, the only 
remaining problem is the convergence of the FET. 


In [Appendix B we formulate and prove the condition which has to be fulfilled by the 


mass matrices in the flavour basis, in order to make EET written for any one-loop amplitude, 
convergent. The result is that the moduli of every eigenvalue of the dimensionless mass 
insertion matrix has to be smaller than one. 


4 Expansion of fermionic amplitudes 

Expanding amplitudes in which flavour violation enters through fermionic mass matrices is 
more complicated, because such matrices do not need to be Hermitian and in general can 
be diagonalized with the use of two different unitary matrices. Nevertheless, as it turns out, 
FET can always apply to this case as well, with minor, but necessary, modifications which 
we discuss below. 

Lets first consider a Lagrangian of A-copies of Dirac fermion free fields. This will have 
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(4.1) 

(4.2) 


the general form, 

4alour = i^A^'ifA-MABi^APL^B) 

in a 4-spinor Dirac notation, which is more suitable for mass eigenstates calculations and 
offers a more compact description in our following discussion. Since Majorana spinors can be 
understood as Dirac spinors with an extra chirality constraint, our discussion applies directly 
to Majorana fermions, as well. 

As is well known, in a chiral theory, a Dirac spinor is in general reducible under flavour 
rotations. The transformation from flavour to mass basis is performed through two different 
unitary matrices, acting independently on its chiral projections, as 

'^LA = UAii^Li , '^RA = VAii^Ri , (4.3) 

which can always bring an arbitrary complex flavour mass matrix M into a real non-negative 
diagonal form, satisfying 

V^MU = m = diag(mi,..., ruAr). (4.4) 

The unitary matrices V and U diagonalize also the (semi) positive-definite Hermitian matrices 
MMt and M^M, through the transformations 

V'fMM'fV = U'^M'fMU = m^. (4.5) 


To streamline the notation, we introduce the unitary matrices lA and W, generalizing our 
transformation rules for chirality projected fermion fields to a reducible Dirac 4-spinor, as 


U = VPl + ^^PR , U = UtpR + \^Pl . 


(4.6) 


In this compact description, eq. (4.3) will result in 

4'a = UAi^l^i, ^A = tpi UiA ■ (4.7) 

The free propagator for the fermion multiplet 'I' in flavour basis is a matrix both in spinor 
and flavour space. Inverting the Dirac operator in eq. (4.2), we find 

i 


A{k) = 


- MPl - MtPj 


R 


= {M^Pl + ^^Pl) 


A:2 - MMt 


+ (MPr + ^Pr) 


- MtM 


(4.8) 


The free propagators in flavour and mass eigenstates basis are related by the same rotations 
as fermion fields. From the identity, 

(0|T{^s(x)^A(y)}|0) = UB^{0\T{^|;i{x)^P^{y)}\0)UiA, (4.9) 

it follows that. 
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where Ai(k) is the fermion propagator in the mass eigenstates basis. Applying the explicit 
expressions of eq. (4.6) for the reducible flavour rotation matrices, to eq. (4.10) and using the 


following algebraic identities: 

d 1 


U' 


U 


k‘2 — 

m 

/j2 _ j^2 


u = 


1 


U = 


A:2 - MMt 
1 


rPL + 


1 


A:2 - MtM 


Pr, 


A;2 - MMt 


Pl + M 


1 


/c2 - MtM 


Pr, 


(4.11) 

(4.12) 


the flavour propagator A.{k) can be also obtained from the mass basis propagator A.{k). 

In order to calculate the amplitude, apart from propagators one needs to consider the 
transformation rules for the fermionic vertices. A general fermionic current in flavour basis, 
can be expressed in the form, 

j'l' = , (4-13) 

where T is an operator acting both in flavour and spinor space, and may also depend on scalar 
or gauge fields. In any QFT model a general fermionic vertex can be decomposed into four 
chirality projected terms as 

f = Prl Pr + PlrPr + PrrPr + Pll Pl 

= PlPrl + PRi'LR + PlPrr PPr^ll ■ (4-14) 

where TiRiRL) a-re scalar- or tensor-type couplings, and FiiiRR) are vector couplings. 

In our compact notation, the transformation rule for vertices can be simply expressed as: 

t = U^TU^, (4.15) 

or explicitly in terms of F components as 

Trl = V Trl Ut , t'RR = Y Trr Vt , 

Tlr = U Tlr vt , tLL = V Tll Ut . (4.16) 


It is important to notice that the transformation rules for propagators and for vertices 
eqs. (4.10) and (4.15), respectively, are different, which reflects the general difference in trans¬ 
formation rules for amputated and non-amputated Green’s functions. 


Let us examine now the general n-point transition amplitude with fermion lines (external 
or internal - our formalism applies to the latter by setting final and initial fermion indices to 
be equal). Lets focus on any chosen fermion line in such an amplitude. The Feynman rule in 
the mass eigenstate basis would have the general form, 


Ai,,~(rAr ...AF)^., . 


(4.17) 


Applying flavour rotation to the external states and using eqs. (4.10) and (4.15)), one can 
get an expression for the fermion line in flavour basis, Mji, written as a sequential product 
of flavour-basis fermion vertices and matrix propagators, 


M 


U^M ~ (U^TU^^ (UAU^ (u^TU^^ ... (uAU^ (u^TU^^ 
f Af ... Af . (4.18) 
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This shows that any fermionic amplitude built of vertices and propagators in the mass basis 
can be formally transformed into the amplitude given in terms of respective quantities in the 


flavour basis. We can now observe that, as shown explicitly in eq. (4.8), matrix denominators 


of the loop integrals in the flavour basis always depend on Hermitian matrices MM^ or 
M+M, and only such combinations would appear as formal arguments of loop functions. As a 
consequence, one can conclude that FET formulated for Hermitian matrices can always apply 
to loop functions appearing in fermionic amplitudes, as well. 

From more practical point of view, our derivation leads to the conclusion that the fermion 
mixing matrices U and V can appear in amplitudes only is some specific combinations, namely 


UBifimi) 

U*M 

= fiM^M)BA, 

Vsi fimf, 


= f{MM^)BA, 

Usi 


= M]j^/(MMt) 

Vsi mif{ni^) 


= MBcfiM^M) 


ft 

‘■CA ’ 


(4.19) 


which can always be expanded using eq. dOI ). 

We should notice that the formal treatment followed in this section can easily generalize 
to the case of more complicated flavour models, where sets of flavour fields belong to distinct 
flavour families. In this case, the propagators and the vertices in the general formulae of 
eqs. (4.17) and (4.18), will carry both internal (flavour) and external (family-group) indices. 
However, only F or F can associate different family groups because propagators, A(/c) or 
A(/c) are block diagonal in family space. Therefore, one can accommodate in this formalism 
amplitudes with complicated flavour structure like, e.g., rare processes in the MSSM with 
fermions on the external lines and sfermions and gauginos circulating in loops. FET formalism 
not only allows one to calculate such diagrams in flavour basis but also in any other “hybrid” 
basis of convenience, e.g., fermions-gauginos in mass and sfermions in flavour basis, or any 
other combination. 


5 Application of FET: neutron EDM in the MSSM 


To illustrate that higher order mass insertion terms can give physically meaningful bounds 
we consider the example of the neutron Electric Dipole Moment (nEDM) in the Minimal 
Supersymmetric Standard Model 

The full nEDM can be expressed as a combination of parton level contributions - EDMs 
of quarks dq, their chromoelectric dipole moments (CDM) Cq and the CDM of gluon Cg. The 
parton moments are defined as respective coefficients in the effective Hamiltonian: 
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ng = (5.1) 


The total neutron EDM depend on its hadronic wave function and can be written as 


Er, 


deddd T Veudu T (^(j]cdC^d T VcuCu) T 


erigAx 

47r 


(5.2) 


12 










where rji and Ax are 0(1) QCD wave-function factors 
scale 
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21 ] and the chiral symmetry breaking 
respectively. Various models can give significantly different values for iji, differing 


even by an overall sign. Thus, eq. (5.2) and the bounds it puts on MSSM parameters should 


be treated as order of magnitude estimates only, since potential cancellations in (5.2) depend 
on these poorly known coefficients. 

The explicit expressions for dq, Cq and Cg are given in ref. [^. In this example, we consider 
only the dominant gluino contribution to the up-quark EDM and CDM. Taken together, their 
contribution to the nEDM can be expressed as 


1 

Ms 


6 




k=l 


Uk‘ 


(5.3) 


where Zjj and are up-squark mixing matrices and physical masses, M 3 is the gluino 


mass (for conventions and the detailed definitions see refs, 
ratios, xq = The function F(x) is the sum of loop contributions 


24,^), and we define the mass 


F{x) = {^r]euCl 2 {x) - (18C'ii(x) -h Ci 2 {x)) 

loTT V Z 


(5.4) 


with C'ii,C'i 2 being certain PV-functions defined as 


Cii(x) 

Cl 2 {x) 


— l-|-3x ^ 

4 ( 1 -x )2 ^ 2(1 -x )3 
X -|- 1 X 

2(1 — x )2 (1 — x )3 


logx, 

■ log X . 


(5.5) 

(5.6) 


Flavour violation in the MSSM squark sector is strongly bounded by numerous experi¬ 
ments and known to be very small, < 0 ( 10 “^), for down squark mass matrices if the diagonal 
elements of those matrices are around the electroweak scale. Therefore, we assume for the 
purposes of this example that the left down soft SUSY breaking squark mass matrix is ap¬ 
proximately diagonal, but not degenerate, of the form 

( 0 0 \ 

0 vr?- F SniL.q 0 . (5.7) 

V 0 0 + / 

In the left up-squark sector the off-diagonal mass terms are then generated by the SU{2) 
relation: 



(rri^(j)LL = K{m?^)LLK\ (5.8) 

where K denotes the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 

Consider now the flavour expansion of eq. (5.3). In the first order it constrains the imag¬ 
inary part of the trilinear up-squark mixing. 



^^^Im(A^i + Y„^Mot/3)FW 


(5.9) 
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where the RHS is now expressed in terms of parameters in flavour basis. In particular, the 
arguments of the first divided difference, are given by diagonal elements of up-squark 
mass matrix, 


^Ul - 


_ 


Mi 




_ 


Mi 


(5.10) 


Following the current experimental bound, \En\ < 2.9 x 10 26 , and bearing in mind 


potential cancellations, eq. (5.9) sets strong bounds on the imaginary phases of // and A]j, of 
the order of 10“^ and 10“^, respectively, for SUSY mass scale of the order of 1 TeV. 

What is interesting, and to our knowledge has not been discussed thus far in the literature, 
is that the experimental bound on nEDM is so strong that it constrains significantly also the 
real parts of up-squark mass insertions, an effect which is easily visible when analyzing higher 
orders in MIA expansion. To avoid lengthy expressions, let us assume that in the up-squark 
sector only the 31 off-diagonal entries do not vanish in the “right-handed” soft mass matrix 
and in the trilinear mixing matrix and that they are purely real. In addition, the 


{rri^)LL is defined by the relation to diagonal down sector in eq. (5.8). Then, a non-vanishing 


contributions to nEDM are generated from higher orders in mass insertions via the mixing 
with the complex CKM matrix elements. Using the FET up to 2nd order one can see that 
the result depends only on the A^, 


^ 2 ) 


D 


V 2 sin 2013 cos 023 

2\/2M| 


(<5m|,i3 - (5m|,i2 sin^ 0 i 2 ) sin 8ckm 








(5.11) 


where 0i2,0i3,023, and 5c km are respectively, the angles and the CP-violating phase in the 
standard CKM matrix parametrization. Note again that a CP-violating observable constraints 
real squark mass parameters through the CKM CP-violating phase. 

It is worth noting that even the 3rd order expansion of FET sets numerically significant 
constraints on the real parts of flavour violating parameters. In particular, the dependence 
on (m|)3i^ parameter, absent at lower orders, is now introduced through. 




D 


V 2 sin 2013 cos 023 


(hm|)i3 — 5mi)i2 ^ 12 ) sin 6 km x 


2V2MI 

X Re(m|)|i^ Re(A^3 + Ytp* cot /3) ,: 


'Ur3 


(5.12) 


Comparing separately expressions given in eqs. (5.11) and (5.12) with the experimental upper 


bound on the neutron EDM, one can obtain order of magnitude estimates on, otherwise weakly 
constrained, 31 and 13 elements of the up-squark trilinear and “right-handed” soft mass terms 
in relation to mass splitting in the down-squark sector. Such bounds are important e.g. for 
analysis of the maximal allowed decay rates of the top quark to lighter MSSM Higgs boson, 
t ^ uh 17 . The numerical results for a typical MSSM parameter set, obtained using the full 
unexpanded mass eigenstates expressions for nEDM and the SUSY_ELAVOR library 27 


31 


are collected in Table They all agree both qualitatively and quantitatively with eqs. (5.11) 
and (5.12) that have been obtained from the FET of eq. (3.11). 


Alternatively, one can satisfy the nEDM bound by choosing flavour violating entries large 


but correlated so that terms in eqs. (5.11) and (5.12) cancel each other to large accuracy. 
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dmois [TeV] 

0.2 

0.4 

0.6 

0.8 

1 

iReA^VMal 

2.43 X 10 "^ 

1.31 X 10 "^ 

1.12 X 10 “^ 

8.94 X 10 “^ 

7.85 X 10 "^ 

|Re(m|)?i^/M|| 

2.47 X 10"2 

1.23 X 10“2 

1.11 X 10"2 

8.64 X 10"3 

7.41 X 10"3 


Table 1: Upper bounds on |ReTy/M 3 I and imposed by current experimental con¬ 
straints from neutron EDM. Displayed values were obtained assuming = 0 for the 2nd row, 

= 0 for the 3rd row and no other sources of the sfermion flavour violation. Other parameters set 
to tan/3 = 4, common SUSY-scale M 3 = 1.1 TeV and a suitable value of was implicitly chosen to 
satisfy the 125 GeV Higgs mass constraint. 
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Figure 3: Range of the KeAfj/Ms — Re(m^)^^^/M| plane allowed by current experimental 
constraints from neutron EDM. MSSM parameters defined as in Table [T| 


It is interesting to observe that CKM-related factors in these equations are identical, so the 
correlation between Af^ and is given only by SUSY parameters (of course once the 

QCD-related factors in eq. (5.4) become better known). Both terms exactly cancel when 


Re[Al^/M 3 ] _ Aff + Ytp* cot j3 




M3 


F[2l I 


( 5 . 13 ) 


Eq. (5.13) suggests a linear shape for the allowed ReA^^ — Re(m^)^^j:j parameter space, 
consistent with the nEDM experimental bound. This result is illustrated clearly in Eig. 
obtained again with the SUSY_FLAV0R code (he., without using the MIA expansion), and as¬ 
suming values of QCD factors implemented there. Again this result follows closely eqs. ( 5.11[ ) 
and (5.12) that have been obtained from the EET. 

Eurthermore, we have successfully applied EET to another observable, namely the rare 
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top decay to light quarks and a Higgs boson, t —)• g/i, in MSSM [^. What one practically 
gains from using FET in complicated theories, like MSSM for example, is to algebraically 
isolate the dominant effects for a given observable without performing tedious diagrammatic 
MIA expansion nor extensive computer scans of a multi-parameters space. 


6 Summary 


In this article we have presented and proved a theorem of matrix analysis, the Flavour Ex¬ 
pansion Theorem (EET), that remarkably translates any transition amplitude, obtained in 
terms of mass eigenstate basis parameters, into its corresponding amplitude in flavour eigen¬ 
state basis, purely algebraically, without the use of standard diagrammatic methods like the 
Mass Insertion Approximation (MIA) method. Following the formulation of this theorem, 
any analytic function of a Hermitian matrix is expanded in terms of its off-diagonal elements 
with coefficients being the multi-variable functions, commonly known as divided differences. 
Natural implementation of such expansion [see eq. (3.11)] comes from the intimate connection 
between the divided differences and the Passarino-Veltman one-loop functions. Apart from 
the formal proof, we have discussed also FET limitations, such as non-analyticity and con¬ 
vergence issues. We have also extended the use of the theorem in case of general transition 
amplitudes involving fermion mass matrices which are not necessarily Hermitian. 


We have argued many times throughout this article, that the algebraic derivation of the 
flavour basis result using EET is substantially easier, shorter and more compact than the 
diagrammatic one. A pedagogical example is given in Section However, we also illustrate 
how FET works with a significant example based on sparticle (gluon-squark) contributions to 
neutron-EDMs. In this example, the use of FET at higher non-trivial orders is capable to set 
fairly strong bounds of order ~ 10“^ on real parts of up-squark mixing matrix elements from 
nEDM measurements by assuming that CP-violation arises only from the CKM-matrix phase. 
This FET result agrees with our exact numerical calculations [see Table and Fig. [^] using 
SUSY_FLAV0R library. To our knowledge these bounds are new in the MSSM flavour physics 
literature and demonstrate the usefulness of the Flavour Expansion Theorem, especially, when 
it applies to complicated models. 
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Appendix A Proof of the Flavour Expansion Theorem 


The FET theorem formulated in Section 3.2 can be proved using mathematical induction and 
the notion of the “fully symmetrized polynomials”. 


A.l Fully symmetrized polynomials. 


The “fully symmetrized” homogeneous polynomials of degree N 32 , can be understood 
through the following equivalent definitions: 

Definition 3. {xi,... ,xm) is the direct sum of all distinct N-degree monomials con¬ 
structed out of the given set of M variables Xi. 

Definition 4. Alternatively, QfJ (xi,..., xm) is defined as 

N 


Q^{xi,...,xm) ^ 


N 


aM — ST 
“^M — V-^1 

ai ,...,aM=0 


aM 

... 


ai-\-...-\-aM—^ 


(A.l) 


Directly from definitions above, one can express the ful ly sy mmetrized polynomial, , 
for any value of M, N. Due to the symmetric nature of eq. (A.l) there exist many equivalent 


representations of this sum. For non trivial M, N, one of these will have the explicit form 

Q JY (xi, . . . , X][/[ ) 

(xi)^ + (xi)^“^(x2 + X3 + . . .) + + . . . + (X2X3) + . . .^ + . . . + Xi(...) 

+ ( 2 : 2 )’^ + (^ 2 )^ ^(X3 + ...) + ... 


(xm) 


N 


(A.2) 


The identities, = 0, = 1, ~ ^ hold trivially by definition. 


Due to eq. (A.l) the factorization property. 


N 


Q^{xi,...,xm) = '^Qk{xi,...,xl)Q^_k{xl+i,---,xm) , 


K=0 


{no sum over L ), 

(A.3) 


holds for any integer L, satisfying 1 < L < M — 1, and for any choice of, {L} and {M — L}, 
respective subsets of M variables, Xj. 


Lemma 1. Fully symmetrized polynomials are M-order divided differences of the gen¬ 
erating functional, thus satisfying Equivalently, the expression, 

QV = Qn+} > (A.4) 


holds for any M > 1. 
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Proof. First, we show the validity of eq. (A.4), for M = 1, namely 


Q% = 


(A.5) 


Applying eq. (A.3), we have 


N 


K=0 


(Xi - X2)Q%{xi,X2) = X] ( (^2) - QK{xi)QN-K+lix2) 


N 


N 


Q]v+i(®i) “ Qn+i{x2) + ^ QN-K+lixi)QK{x2) — ^ Qk{xi)Qn-K+i{x2) 


K=1 


K=1 


— Qn+Axi) — Qn+Ax2) 


(A.6) 


which is equivalent to eq. (A.5). Now it is straightforward to verify eq. (A.4) for M > 1, as 
well. Denoting y = {xs,..., xm+i}, we have 

,_, N 

fO 


(xi - X 2 )Q^^^{xi,X 2 ,y) ^ (xi - X 2 ) QK{xi,X 2 )Q^_K{y) 

K=0 
N 

E {QK+li^A-Q]<+l{x2))QYKiy) 


K=0 

N+1 


E {QKixi) - QA{x2)'^QN-YAy) + (< 30 ( 3 ^ 1 ) - Qlix2))QN+liy) 


K=1 

iM 


= Qjv+1 {xi,y)- Qn+1 {x2 , v ), 

and therefore finishing the proof of the lemma. 


(A.7) 

□ 


A.2 Flavour Expansion Theorem: the proof 

Proof. Due to Def.the Hermitian matrix function /(A), can be expressed as a power series, 

00 

/(A) = E^-A-. (A.8) 

m=0 

One can apply the matrix decomposition A = Aq + A to the above series (convergent by 
assumption) and rearrange terms collecting together the same powers of A. Assuming that 
the resulting summation remains convergent, we have 

/(A) = Co I + ciAo + C2 Aq + C 3 Aq + ... : Fq 

+ ciA + C2 (AAq) + C 3 (AAq) + C 4 (AAq) + ... : Fi 

+ C2 A^ + C3 (A^Ao) + C 4 (A^Aq) + ... : F 2 


+ + cm+i(A^Aq) + ... + c„(A^Aq + ... : Fm 

+ cm+iA^+i + ... + + ... : Fm +1 

(A.9) 
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where we have defined 


(A”^AS) = ^ A”^AS , (A.IO) 

P—distinct 

for all distinct permutations of the set {A ..., Aq ...} of (m + n) objects. 

The matrix element (/|/(A)| J) will be given by the sum 

OO 

(/(^))'^= + H'' + ■■■ , (A.ll) 

N=0 

where, by direct calculation, the above terms are given by (summation over repeated internal 
indices Aj, is considered - also if they appear more than twice), 

, (A. 12 a) 

fI^ = Aij[ci + C 2 [Aq + Aq] + C 3 [(Aq)^ + (Ag)^ + AqAq] + ... ^ , (A.12b) 

= AjK^Ak^j{^2 + C 3 Aq + Ag + A^^ 

+ Ci[{Aif + {Aif + {A^^f + AiAi + AiA^^+AiA^^^^ + ...) 

OO 

= Aik^AkiJ ^ C 2 +Ar( 5 Ar(Ao, Aq, A^i) , (A.12c) 

Ar =0 

OO 

Fm = Ajk^Ak^K2 ■ ■ ■ CM+jvQnA q, Ao,A^\ ■ ■., A^^~^) , (A.12d) 

Ar =0 

To prove the theorem, we need to show 

OO 

CM+NQT\^i^ • • • > (4, Ai, A^\..., A^--) , (A.13) 

N=0 

for all M > 0. This can be realized using mathematical induction. For M = 0, we obtain 
trivially 

OO OO 

E = E = nA) ^ ■ ( a .14) 

N=0 N=0 

Now, let 

OO 

^ cm+nQn^^ = , (A.15) 

Ar =0 

holds for some M > 0 and for any set of M + 1 arguments. Then, we need to show that 

OO 

CM+N+iQ^^\Ai, Ai, A^\..., A^-) = Ai, A^S ..., A^-) , (A.16) 

Ar =0 
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which, by Def. [^of divided differences in eq. (3.7), is equivalent to showing 

oo 

{Aq — Aq ) ^ cm+n+1 (^ 0 ) ^ 0 ) ^0 ^) • ■ ■) ^0 ^) 

N=0 

_ __ oo 

^ £ cm+n+1 • • 


iV=0 


oo 

CM+N{Q^^HAi,A^\.. .) - Q^+\AiA^\.. .)) 

iV=l 

+ CM{Q^^\Ai • • •)) 

oo 

= ^ cm+n (Aq, Aq^ , ...) — Q^~^^{Aq,Aq^ ,.. .)^ 

N=0 

and hence the theorem is proved. 


(A.17) 

□ 


Appendix B Convergence criterion for FET expansion of the 
one-loop functions 


It is well known that, any one-loop amplitude can be expressed as a linear combination of 
“master” PV-integrals with trivial i.e., equal to 1, integrand numerator. Thus, it is sufficient 
to find a convergence criterion for the FET expansion only for master integrals. Below we 
formulate such a criterion for the most often considered case of loop functions with vanish¬ 
ing external momenta. The same criterion can be applied to coefficients of the expansion 
of one-loop integrals in terms of external momenta (assuming that they are far from thresh¬ 
olds and momentum expansion can be performed) - such coefficients can be also reduced to 
combinations of master integrals with vanishing momenta. 

For vanishing external momenta master integrals can be expressed as 


PVt\ml 




—z(47r)^ 


d^p 1 



udu 

n;=i(«+ruj) 


(B.l) 


where we assume n > 3 to avoid divergent integrals - considering the estimates for finite ones 
is sufficient to establish the convergence criterion for FET expansion as it depends only on 
behaviour of higher order terms. 


Eq. (B.l) leads immediately to inequality 




^ -A 
m 


n-hl 


... ,ml) 


(B.2) 


Applying this inequality, iteratively for higher order terms, to majorize the RHS of eq. (3.14) 
(in what follows we do not write explicitly any PV-function arguments apart from the ones 
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used in the expansion), we get 


Pvt\A) 


J IJ 


< 


< 


6jjPvi^\Ai) + Pvi^+^\Ai,Ai)Aij 


+ 


PVq {Aq,Aq,Aq)AjkAkj 


+ ... 


6ijPV^^\Ai) + PV^^+^\Ai,Ai) Ajj 


+ 


< 


'Aij 


Aik 


pvriK) 


Pvt\Ai) 


5IJ + 




+ 


Akj 
Aik 


+ ... 

Akj 






+ .. 




Aij 

- + - 

Aik 


(B.3) 


Akj 


»V\/^ 


+ ... 


where we assume that all indices apart from I, J are implicitly summed in the range 1... N. 
Let us now define the symmetric matrix Q with elements being the absolute values of dimen¬ 
sionless quantities commonly referred in the literature as “mass insertions” (diagonal elements 
of Q vanish by definition of A matrix) 



Then eq. (B.4) can be expressed as 





< 


pvt\4) 




(Q + Q2 + .. 



(B.4) 


(B.5) 


The expression in the inner parenthesis of the RHS of eq. (B.5) is a geometric series. According 
to the definition of a function of Hermitian matrix given in Section 3.1, this series is convergent 
if it converges also for any of Q eigenvalues, hence their absolute values must be all smaller 
than 1. This can be expressed formally, as 


sup je^Qel = sup je^Dgel < 1, (B.6) 

||e||=l ||e||=l 

where e denotes any real unit vector, and Dg is the diagonal matrix of eigenvalues. Obviously, 
this is a sufficient but not necessary condition for the convergence of the expansion. 

Finally we should note that vanishing diagonal elements can not pose a threat for the con¬ 
vergence of the FET expansion in physical theories. This is because all Hermitian (squared) 
mass matrices are semi-positive definite matrices, and for such matrices if Ag = 0, then 
necessarily also Aik = Aki = 0 for all K. Thus all potentially divergent terms vanish. 
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